UNCLASSIFIED _ 

SECURITY  CLASSIFICATION  OF  THIS  PAGE 


REPORT  DOCUMENTATION  PAGE 


la.  REPORT  SECURITY  CLASSIFICATION 

Unclassified 


lb.  RESTRICTIVE  MARKINGS 

None 


2a.  SECURITY  CLASSIFICATION  AUTHORITY 


3.  DISTRIBUTION  AVAILABILITY  OF  REPORT 


2b.  DECLASSIFICATION  DOWNGRADING  SCHEDULE 


Approved  for  public  release;  distribution  is 
unlimited. 


4.  PERFORMING  ORGANIZATION  REPORT  NUMBER(S) 

NORDA  Report  1 34 


5.  MONITORING  ORGANIZATION  REPORT  NUMBER(S) 

NORDA  Report  134 


6.  NAME  OF  PERFORMING  ORGANIZATION 

Naval  Ocean  Research  and  Development  Activity 


7a  NAME  OF  MONITORING  ORGANIZATION 


Naval  Ocean  Research  and  Development  Activity 


6c.  ADDRESS  (City,  Stale,  and  ZIP  Code) 

Ocean  Science  Directorate 
NSTL,  Mississippi  39529-5004 


7b  ADDRESS  (City,  State,  and  ZIP  Code) 

Ocean  Science  Directorate 
NSTL,  Mississippi  39529-5004 


8a.  NAME  OF  FUNDING 'SPONSORING  ORGANIZATION 

Defense  Mapping  Agency 


8b  OFFICE  SYMBOL 
(If  applicable) 


9  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 


8c  ADDRESS  (City,  State,  and  ZIP  Code) 


10  SOURCE  OF  FUNDING  NOS. 


PROGRAM 

PROJECT 

TASK 

WORK  UNIT 

Washington,  D.C. 

ELEMENT  NO 

NO. 

NO 

NO 

1 1  TITLE  (Include  Security  Classification) 

Sounding  Vessel  Position  from  Adjustment  by  Variation  of  Parameters 


1  2  PERSONAL  AUTHOR(S) 

G.  Mapp,  H.  J.  Byrnes,  and  M.  Heinzen* 


13a  TYPE  OF  REPORT 

13b.  TIME  COVERED 

14  DATE  OF  REPORT  (Yr.,  Mo.,  Day) 

15  PAGE  COUNT 

Final 

From  To 

December  1985 

15 

16.  SUPPLEMENTARY  NOTATION 

*  Naval  Oceanographic  Office 


1 7  COSATI  CODES 

FIELD 

GROUP 

SUB  GR. 

18  SUBJECT  TERMS  (Continue  on  reverese  if  necessary  and  identify  by  block  number) 

coastal  survey,  navigation,  least  squares,  adjustment 


1  9.  ABSTRACT  (Continue  on  reverse  it  necessary  and  identify  by  block  number) 

For  hydrographic  surveys  conducted  by  the  Naval  Oceanographic  Office,  the  position  of  the  sounding 
vessel  is  determined  by  applying  the  method  of  adjustment  by  variation  of  parameters.  Three  types  of 
navigational  aids  are  used;  ranging  (ARGO,  RAYDIST,  LORAN-C),  azimuthal  (Coast  Artillary  Azimuth  In¬ 
strument),  and  hyperbolic  (ARGO,  DECCA,  HI-FIX6,  LORAN-C).  Given  data  from  any  combination  of  at 
least  two  navigational  aids,  a  fix  may  be  obtained  using  an  iterative  method,  which  applies  successive 
adjustments  to  an  approximate  location  and  forces  it  to  converge  to  the  most  probable  position.  The 
magnitude  and  direction  of  each  adjustment  is  determined  from  a  least  squares  solution  that  minimizes 
the  residual  differences  between  actual  navigational  observations  and  imaginary  observations  calculated 
as  if  the  ship  were  at  the  approximate  location. 


20.  DISTRIBUTION/AVAILABILITY  OF  ABSTRACT 
UNCLASSIFIED/UNLIMITED  □  SAME  AS  RPT. 


DTIC  USERS  □ 


21  ABSTRACT  SECURITY  CLASSIFICATION 

Unclassified 


22a  NAME  OF  RESPONSIBLE  INDIVIDUAL 

H.  J.  Byrnes 


22b.  TELEPHONE  NUMBER  (Include  Area  Code) 

(601)  688-5257 


22c.  OFFICE  SYMBOL 

Code  352 


DD  FORM  1473,  83  APR 


EDITION  OF  1  JAN  73  IS  OBSOLETE. 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PAGE 


LIBRARY 

RESEARCH  REPORTS  DIVISION 
NAVAL  POSTGRADUATE  SCHOOL 
MONTEREY,  CALIFORNIA  93940 

Naval  Ocean  Research  and  Development  Activity* 

^December  1985  Report  *134  . 


Sounding  Vessel  Position  from 
Adjustment  by  Variation  of  Parameters . 


George  Mapp 
H.  J.  Byrnes 

Mapping,  Charting,  and  Geodesy  Division 
Ocean  Science  Directorate 

Morris  Heinzen 

Naval  Oceanographic  Office 


Sponsored  by: 

Defense  Mapping  Agency 
Washington,  D.C. 


Approved  for  public  release;  distribution  is  unlimited.  Naval  Ocean  Research  and  Development  Activity,  NSTL,  Mississippi  39529-5004. 


Foreword 


This  report  is  the  first  of  a  series  that  will  provide  documentation  for  soft 
ware  being  developed  specifically  for  Naval  Oceanographic  Office  hydrographic 
survey  applications.  New  hardware,  software,  and  procedures  that  maximize 
automation  to  improve  efficiency  while  maintaining  high  quality  of 
hydrographic  products  are  being  treated  collectively  as  a  Hydrographic  In¬ 
formation  Handling  system.  Future  reports  will  cover  various  aspects  of  the 
total  processing  system. 


R.  P.  Onorati,  Captain,  USN 
Commanding  Officer,  NORDA 


Executive  summary 


For  hydrographic  surveys  conducted  by  the  Naval  Oceanographic  Office, 
the  position  of  the  sounding  vessel  is  determined  by  applying  the  method 
of  adjustment  by  variation  of  parameters.  Three  types  of  navigational  aids 
are  used:  ranging,  azimuthal,  and  hyperbolic.  Given  data  from  any  combina¬ 
tion  of  at  least  two  navigational  aids,  a  fix  may  be  obtained  using  an  iterative 
method,  which  applies  successive  adjustments  to  an  approximate  location  and 
forces  it  to  converge  to  the  most  probable  position.  The  magnitude  and  direc¬ 
tion  of  each  adjustment  is  determined  from  a  least  squares  solution  that 
minimizes  the  residual  differences  between  actual  navigational  observations 
and  imaginary  observations  calculated  as  if  the  ship  were  at  the  approximate 
location. 
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Sounding  vessel  position  from  adjustment 
by  variation  of  parameters 


•  I.  Introduction 

The  Naval  Oceanographic  Office  conducts  bathymetric 
surveys  as  part  of  its  continual  effort  to  chart  the  world 
ocean.  For  this  work,  it  is  essential  to  determine  sound¬ 
ing  positions  with  a  precision  compatible  with  the  product. 
^  Coastal  survey  work  is  done  using  three  types  of  naviga¬ 

tional  aids  (navaids):  range,  azimuth,  and  hyperbolic  (Figs. 
1-3).  Here,  a  hyperbolic  navaid,  wThich  is  really  a  pair 
of  ranging  stations,  is  considered  to  be  a  single  navaid. 
From  each  navaid  is  obtained  a  line  of  position;  the  in 
tersection  of  lines  of  position  determine  the  ship  position. 
®  If  there  are  more  than  two  lines  of  position,  the  intersec¬ 

tion  generally  is  not  a  point.  In  this  event,  the  most  prob¬ 
able  position  is  found  by  weighting  each  navaid  according 
to  its  precision  and  selecting  the  position  accordingly. 

The  method  presented  here  differs  from  other  published 
^  work  in  that  combinations  of  different  navaid  types  are 

used  to  obtain  a  fix. 


II.  Adjustment  by  variation 
of  parameters 

Given  data  from  any  combination  of  two  or  more 
navaids,  the  most  probable  position  may  be  calculated  by 
an  iterative  method.  First,  the  approximate  location  is 
specified.  For  this  position  a  set  of  imaginary  observations 
is  calculated.  A  new  approximate  position,  which  has  im 
aginary  observations  closer  to  the  actual  observations,  is 
then  calculated.  The  process  of  adjusting  an  approximate 
position  is  repeated  until  the  magnitude  of  adjustment 
becomes  sufficiently  small,  at  which  time  the  approximate 
position  is  taken  as  the  actual  position. 

The  terminology  encountered  in  the  literature  is 
somewhat  confusing.  It  is  important  to  distinguish  between 
an  actual  observation,  which  is  a  measurement  from  a 
navaid,  and  an  imaginary  observation,  which  is  a  value 
that  one  would  expect  to  measure  at  a  given  location.  Im¬ 
aginary  observations  are  often  referred  to  as  “computed 
observations’ ’  because  the  source  of  the  “observations” 
is  a  calculation  that  uses  geodetic  formulas  to  compute 
a  given  range,  range  difference,  or  azimuth. 


Figure  1  Azimuth  navaid  Transit  or  theodelite  is  used  for  this 
measurement.  Azimuth  is  the  sum  of  the  offset  and  observa¬ 
tion  angles;  this  sum  yields  the  bearing  of  the  ship  with  respect 
to  the  center  station .  The  line  of  position  is  a  straight  line. 

A.  The  general  model  of  variation 
of  parameters 

The  derivation  of  the  general  model  for  variation  of 
parameters  follows.  A  more  thorough  discussion  of  this 
model  may  be  found  in  Ewing  and  Mitchell  (1970)  and 
Mikhail  and  Grade  (1981).  Initially  there  are  two  data 
items:  the  approximate  position  (jr,  y)  and  the  actual  navaid 
observations  (/.).  The  first  step  is  to  calculate  the  im¬ 
aginary  observations  (I A  using  geodetic  formulas  (e.g.,  the 


Figure  2.  Range  navaid.  Ship  emits  pulse  that  is  simultaneously 
received  and  retransmitted  by  a  transponder.  Round  trip 
distance  is  proportional  to  the  elapsed  time  between  transmis¬ 
sion  and  reception  by  the  ship.  The  line  of  position  is  a  circle. 
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Figure  3.  Hyperbolic  navaid.  Master  sends  out  pulse ,  which 
is  received  by  the  ship  and  slave.  After  waiting  for  a  specified 
period  (the  delay),  the  slave  emits  a  new  pulse.  The  ship 
measures  the  difference  in  time  of  receipt  of  master  and  slave 
pulses.  Knowing  the  baseline  distance  (this  is  published 
information)  the  navigator  can  calculate  the  difference  in 
distance  of  the  ship  from  the  master  and  slave.  The  line  of 
position  is  a  hyperbola ,  the  locus  of  points  having  the  same 
difference  in  distance  between  the  master  and  slave.  Actual 
ly  the  above  is  an  oversimplification  of  the  method  used. 
The  signals  of  the  master  and  slave ,  rather  than  being  pulses, 
are  continuous  sine  waves.  The  ship  detects  the  phase 
difference  between  the  master  and  slave ;  lane  count  is 
kept  by  the  instrument  and  must  be  initialized  at  a  known 
location. 


The  subscripts  /  and  j  are  navaid  number  and  iteration 
number,  respectively. 

The  goal  is  to  calculate  the  appropriate  adjustment  that 
will  “move’'  the  approximate  location  in  such  a  way  that 
the  difference  between  the  actual  and  imaginary  observa¬ 
tion  is  reduced.  After  adjustment,  a  new  set  of  imaginary 
observations  (/  .+  })  may  be  calculated.  The  quantity  /  ;+ 1 
may  be  expressed  in  terms  of  the  previous  set  of  imaginary 
observations  (L.)  and  a  small  change  in  these  values 


Vi  '  '«  +  d'» 


(2) 


or  in  terms  of  the  actual  observations  (/.)  and  the  devia¬ 
tion  from  actual  observations  (v.). 


i,j  + 1  i  i 


(3) 


The  quantity  v.  is  referred  to  as  the  residual  difference. 
Similarly,  the  difference  between  actual  and  imaginary 
observations  is  the  misclosure  (/■). 


A  =  h  -  hi 

Rearrange  Equation  (3)  to  yield 


(4) 


v.  =  /.  -  I  . 

t  t  :,j+ 1 


(5) 


Now,  substitute  values  of  /  -+J  and  /■  from  Equations  (2) 
and  (4) 


vi  -  ft  +  V  -  ft,  + 

V.  m  /.  -  dl.  . 

I  J  l  l] 


Expand  the  dltj  term: 

dF.  (x,  y)  dF-  (x,  y) 

dl.  =  — — —  dx  +  — - —  dy 

’’  d  x  dy 


(6) 


(7) 


Substitute  this  expression  for  dL.  into  Equation  (6). 
dF  dF. 

vi  =  fi  -  -f  dx  -  -f  dy  (8) 

ox  dy 


This  is  the  observation  equation,  and  there  is  one  for  each 

Sodano  inverse  method,  Campbell,  1964).  Imaginary  obser-  navaid.  The  set  of  observation  equations  may  be  expressed  # 

vations  depend  on  the  approximate  location.  concisely  in  matrix  form  as  follows. 


hi  -  Fi  fc  y) 


(!)  v  =  f  -  B  A  , 


(9) 


2 


where 


II 

<1 

V  = 

vl 

V2 

it 

h 

f2 

V 

n 

fn 

B  = 


Since  4>  is  a  scalar,  each  term  on  the  right  side  of  Equation 
(13)  is  also  a  scalar.  The  transpose  of  a  scalar  is  equal  to 
itself.  Therefore,  the  second  term  is  equal  to  the  third. 

A  WWf  =  (biB'Wf)1  =  fWBA  (14) 

So  Equation  (13)  reduces  to 
0  =  flWf  -  2flWBA  +  A'B'WBA  .  (13) 

To  find  the  appropriate  adjustment  A  that  yields  the 

minimum  value  of  0,  set  the  partial  derivative  —  equal 

3A 

to  zero  and  solve  for  A.  This  is  done  in  Equations 
(16)-(20).  See  Mikhail  and  Gracie  (1981,  pp.  73  and  322) 
for  calculation  of  the  following  derivative. 


The  next  task  is  to  find  an  expression  for  the  sum  of 
squares  of  residuals  (0).  For  observations  of  equal  precision, 


0  =  £v  > 


(10) 


but  this  is  rarely  the  case.  More  generally  observations 
are  of  unequal  precision  and  are  weighted  accordingly. 


^  =  2f‘WB  +  2A!(BlWBj  =  0 
dA 


Rearrange  and  divide  by  2. 
A  ‘(B‘WB)  =  f‘WB 


(16) 


(17) 


</>  =  v> 

In  matrix  form,  this  equation  is 


(11) 


Again,  since  each  side  of  the  equation  is  a  scalar,  each 
side  may  be  transposed. 


<t>  ■ 

=  v‘ 

W  v 

(12)  (BWB)‘A  =  B‘W‘f 

where 

B‘WBA  =  B‘Wf 

W, 

0  .  . 

.  0 

Now,  let  N  =  B'WB 

(18) 

W  = 

0 

w2  . 

.  0 

and  /  =  B'Wf  - 

(19) 

0 

0  .  . 

.  w 

n 

then,  N A  =  t 

and  A  =  N~]t  . 

(20) 

Now,  substitute  the  value  of  v  from  Equation  (9). 

<t>  =  (f  -  BA)1  W(f  -  BA) 

=  (f  -  A *B*)  W(f  -  BA) 

=  (fW  -  A tBtW)  (f  -  BA) 

=  fWf  -  A  !BWf  -  flWBA 
+  AtBtWBA 


(13) 


We  have  achieved  our  goal;  Equation  (20)  yields  the 
adjustment  to  the  approximate  position.  The  process  of 
adjustment  is  repeated  until  the  magnitude  of  adjustment 
is  sufficiently  small. 

To  solve  Equation  (20),  one  must  be  able  to  evaluate 
the  j B,  f  and  W  matrices.  Expressions  for  these  quantities 
are  derived  in  the  next  section  for  ranging,  hyperbolic, 
and  azimuthal  navaids. 
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B.  Derivation  of  formulas  for  the 
B,  f  and  W  Matrices 

These  derivations  closely  follow  the  work  of  Heinzen 
(1977). 

1.  B  matrix 


We  start  with  \)  and  derive  expressions  for 
dF  dF 

—  and  — .  Here,  0  and  X  are  latitude  and  longitude, 

30  d\ 

respectively,  and  replace y  and  x  previously  used.  We  will 
derive  expressions  for  range,  hyperbolic,  and  azimuth 
navaids. 


a .  Range 

Here  F(<p ,  \)  =  s,  the  distance  from  navaid  to  ship. 
For  the  first  part  of  this  derivation,  the  shape  of  the  earth 
is  approximated  by  a  sphere.  The  problem  is  solved  by 
constructing  the  spherical  triangle  of  Figure  4.  (See  Ewing 
and  Mitchell,  1970,  p.  295-296  for  geometrical  formulas 
used  here.)  The  subscript  0  indicates  ship;  1  indicates 
navaid.  We  will  switch  back  and  forth  between  the  tw7o 
sets  of  notation  used  in  Figure  4. 

Start  with  the  law  of  cosines. 


cos  a  =  cos  b  cos  c  +  sin  b  sin  c  cos  A 
cos  sQ]  =  cos(90  -  (f>j)  cos(90- <f>0)  +  sin(90- <f)}) 
sin(90  -  4>0)  cos(\1  -  \) 
cos  sQ1  -  sin  0;  sin  <f>0  +  cos  07  cos  0 

cos(\  -  \Q)  (21) 

Note  that  sQJ  is  an  angle  here.  Next,  differentiate 
—  sin  sQ}  dsQ1  =  cos  07  sin  0,;  df>1 
+  sin  0;  cos  0  dcf) 

-  sin  4>}  cos  0  cos(\1  -  \Q)  d4>} 

-  cos  0 j  sin  00  cos(\}  -  \Q)  d(j){ 

-  cos  <f>j  cos  0O  sin(\}  -  \Q)  d(\}  -  \Q)  (22) 

In  the  above  equation  (0;,  X;)  is  the  location  of  the 
transponder  (or  shore  transmitter),  which  is  invariant. 
Therefore, 

d<t>}  =  d\}  =  0 

-sin  sQ1  dsQ]  =  j'sin  0;  cos  <t>0 

-  cos  (f>}  sin  00  cos(\j  -  \Q)J  dcj>0 

+  cos  <t>j  cos  0O  sin(\]  -  \QJ  d\Q  .  (23) 


The  first  term  on  the  right,  ignoring  the  dcj)  factor,  is 

sin  0;  cos  00  -  cos  0;  sin  (f)Q  cos(\}  -  \Q)  .  (24) 

Using  the  following  identity 

cos  b  sin  c  —  sin  b  cos  c  cos  A  =  sin  a  cos  B  , 

Equation  (24)  can  be  written  as 

cos(90  -  <t>})  sin  (90  -  <j)Q)  -  sin(90  -  <j>}) 

cos(90  -  0O9  cos(\}  -  \Q)  =  sin  sQ]  cos  aQ].  (25) 

The  second  term  on  the  right  of  Equation  (23),  again  ig¬ 
noring  the  d\f  factor,  is 

cos  0;  cos  00  sin(\j  -  \Q) 

=  sin(90  -  0;;  sin(90  -  <f>0)  sin(\  -  \Q) 

=  sin  b  sin  c  sin  A .  (26) 

From  the  law  of  sines  we  know 

sin  a  sin  c  ,  .  .  ^ 

-  =  -  and  sin  c  sin  A  =  sin  a  sin  C  . 

sin  A  sin  C 


A  (NORTH  POLE) 


B  (SHIP) 


center  of  earth 

a 

“  S01 

North  Pole 

b 

=  90-*, 

Ship 

c 

=  90-*o 

Navaid 

A 

=  A1  -  A, 

latitude,  longitude  of  ship 

B 

=  °01 

latitude  longitude  of  navaid 

C 

=  a10 

Figure  4 .  Spherical  triangle  used^  in  derivation  of  B  matrix 
expressions.  Note  that  arcs  AB  and  AC  are  meridians: 
therefore ,  b  =  90- 0;,  c  =  90- 0O,  and  A  -  X7-X0. 


4 


Therefore,  Equation  (26)  may  be  written  as 

#  sin  b  sin  c  sin  A  =  sin  b  sin  a  sin  C 

-  sin(90  -  <t>2 )  sin  sQ]  sin  a}0 
=  cos  <j>}  sin  s01  sin  ot}()  (27) 

Substituting  Equations  (25)  and  (27)  into  Equation  (23) 
®  yields 


—  sin  sQ}  ds0J  =  sin  sQ}  cos  a. }]  d<\>{) 

+  cos  <t>}  sin  sQ]  sin  a]0  d\ 

dsQ]  =  -cos  aQ1  d<j>0  -  cos  <t>}  sin  a}0  d\Q  (28) 


To  convert  the  angle  sQ1  into  distance  (arc  length),  multi¬ 
ply  by  radius  of  the  earth.  We  drop  the  subscripts  from 
s  here,  and  shift  from  the  sphere  to  the  spheroid  as  our 
model.  Accordingly,  the  radius  in  the  north-south  dimen- 

#  sion  is  the  mean  radius  of  curvature  (pQ)  in  the  plane 
of  the  meridian;  in  the  east-west  dimension  it  is  the  radius 
of  curvature  in  the  prime  vertical  (j^). 

ds  =  —  p  g  cos  otQ1  d<j>0  -  vQ  cos  <t>}  sin  aw  d\  (29) 

#  Therefore, 


ds 

~  ,  P  o  cos  aoi 

d<t>n 


(30a) 


ds  . 

—  "  -  vo  C0S<t>l  Stn(X10  • 


(30b) 


The  line  connecting  the  ship  and  the  navaid  is  a  geodesic 
line  (Fig.  5).  Therefore,  the  azimuth  a  }]  (navaid  observed 
from  ship)  and  back  azimuth  aJ0  (ship  observed  from 
navaid)  are  related  as  follows. 

cos  <t>0  sin  ocQ]  =  cos  <f>}(-sin  a]Q) 

Now,  Equation  (30b)  may  be  expressed  as 


d s  .  . 

TT-  =  v0  cos  <t>0  Sln  a01 

o\0 


(3r 


Finally,  Equations  (30a)  and  (31)  are  rewritten  for  azimuth 
a0I  to  be  expressed  as  bearing  A(J]  clockwise  from  N. 


dr 

d<A0 


-  P0  cos  aoi 
Po  cos  A01 


~  P0  cos(-A01)  = 


dr 

d\0 


V0  cos  (f>0  sin  a0] 


v0cos  <f)0sin  A01  (32b) 


where  r  =  distance  between  navaid  and  ship. 

p0  =  radius  of  curvature  in  the  plane  of  the 
meridian, 

vQ  =  radius  of  curvature  in  the  plane  of  the  prime 
vertical, 

AQ]  =  azimuth;  bearing  of  navaid  as  observed  from 
ship, 

00  =  latitude  of  ship,  and 
X(  =  longitude  of  ship. 

b.  Hyperbolic 

Here  F(<j),  \)  =  d  =  Rm  -  Rr,  the  difference  in  range 
between  master  and  slave.  The  derivation  is  analogous 
to  the  previous  one,  yielding 


J-  =  p0(cos  A 

OCPo 


01m 


cos  Aou) 


(33a) 


dd 

d\0 


V0  C0S  <t>0(Stn  A01m  ~  Sin  AoJ 


(33b) 


where  the  subscripts  m  and  s  stand  for  master  and  slave. 

c.  Azimuth 

Here,  F((f>,  \)  =  A  =  bearing  of  ship  as  observed  from 
navaid  station.  This  derivation  is  in  plane  coordinates, 
which  suffices  because  azimuth  measurement  may  be  done 


Figure  5.  The  geodesic  line  has  the  following  property: 
cos  <fc0  sin  aQ  =  cos  <f>j  sin  a}  =  cos  </>?  sin  et7, 

(After  Ewing  and  Mitchell,  1970 ,  Fig.  4-6.) 


(32a) 
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only  within  a  short  distance  of  the  ship  (the  line  of  sight). 
Refer  to  Figure  6. 


Next,  substitute  Equations  (35a)  and  (35b)  into  Equation 
(34)  and  calculate  the  following  partials. 


s  —  (x2  +  y2)l/2  =  distance  from  navaid  to  ship 


A  =  arc  tan  ^ 

dA  =  — — - where  U  =  —  and  dU  —  —  — —  dy 

i  +  u2  y  y  y2 


1  ,dx  £,  \  _  ydx  xdy 


l+x2/y2\y  y2 


2,2  2,2 
y  +  x~  y  +  x 


=  (ydx  -  xdy) 
s2 


=  -  (cos  A  dx  -  A 


(34) 


Now,  convert  from  plane  to  spherical  geometry. 


d A  ^  sin  A 

90  Po  s 

(36a) 

— —  =  —  coscf)  cos  A 
d\  s  0 

(36b) 

2.  Weight  matrix 

Expressions  for  computing  weight  are  empirically  de¬ 
rived.  The  basic  principles  are  simple.  Azimuth  weights 
are  constants  and  depend  solely  on  the  precision  of  the 
instrument.  Range  weight  depends  on  instrument  calibra¬ 
tion  and  the  distance  to  the  transponder  (precision  decreases 
with  distance  from  the  transponder). 

Generally,  weight  is  derived  from  the  following  equation. 


_ standard  variance _  # 

variance  of  this  type  of  navaid 


dx  =  east-west  distance 

=  r  cos  4>0  d\  =  v(j  cos  4>0  d\  (35a) 

dy  —  north-south  distance 

=  r  d<f>  =  pQ  d(j)  (35b) 


An  appropriate  value  for  standard  variance  has  been  deter 
mined  to  be  16  m2,  which  corresponds  to  a  standard 
deviation  of  4  m.  As  mentioned,  navaid  variance  for  an 
azimuth  instrument  is  constant  (e.g.,  a  =  0.01°  for  the 
Coast  Artillery  Azimuth  Instrument).  For  ranging  in¬ 
struments,  variance  is  determined  from  this  expression. 


where  o(  -  standard  deviation  of  the  instrument, 
s  **  distance  from  navaid  to  ship, 
m  =  meters. 


The  first  term  is  variance  attributable  to  the  precision  of 
the  instrument.  The  second  term  accounts  for  that  por¬ 
tion  of  variance  that  is  proportional  to  the  square  of  the 
distance. 

Hyperbolic  weight  is  treated  in  a  similar  fashion  to  range; 
the  difference  is  that  the  calibration  errors  of  the  two  com- 
ponent  ranges  cancel  each  other. 


o 2  =  2a2  +  ( - ^  m2 

0  \  10  km  J 

where  oQ  =  (0.02  fx sec)  (299.67  mlfisec) 
d  =  range  difference 
km  =  kilometers. 


Figure  6.  Plane  triangle  used  for  derivation  of  B  matrix  ex¬ 
pression  for  the  case  of  an  azimuth  navaid. 
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Note  that  the  units  for  weights  are  different  for  azimuth 

2  m2 

( — — — \  and  range  ( — A  .  This  is  appropriate  because 
\radian2 )  \m2  ) 

weights  are  multiplied  by  residuals,  yielding  the  result  <f> 
in  meters2. 

3.  Misclosure 

Misclosure  is  the  difference  between  actual  and  im¬ 
aginary  observations,  and  is  calculated  from  the  follow¬ 
ing  expressions. 

a.  Azimuth 

f  =  aim  -  Uobs  =  Aim  ~  (A offset  +  A  obs)  ' 

where  the  subscript  im  =  imaginary  and  obs  =  observed; 
A  =  bearing  from  north. 

b.  Range 

NX 

—  ^ im  ^  obs  ^  im  ^  * 

where  N  =  number  of  lanes  and  X  =  lane  width. 

c.  Hyperbolic  (see  Fig.  3g). 

f  =  dim  ~  dobs 

where  dQb$  -  C  (Atime  -  delay) 

C  =  speed  of  light 

d  =  imaginary  observation  of  range  difference 

R  -  (Rr  +  RJ 

Rb  =  baseline  distance 

Baseline  distance  is  included  in  these  expressions  of  range 
difference;  it  cancels  out  when  d^  is  subtracted  from 

d.  . 

tm 

III.  Discussion 

This  method  can  be  conveniently  implemented  using 
a  computer  program.  An  algorithm  is  described  in  Ap¬ 
pendix  I. 


The  user  must  supply  an  approximate  position  for  each 
fix.  This  approximate  position  is  subject  to  two  constraints: 
it  should  be  close  to  the  actual  position  (within  approx¬ 
imately  100  km);  and  if  either  a  hyperbolic  navaid  or  a 
pair  of  ranging  navaids  is  used,  the  approximate  point  must 
be  on  the  correct  side  of  the  baseline  connecting  the  pair 
of  stations. 

For  a  series  of  fixes  from  a  small  survey  area,  an  initial 
position  may  be  specified,  and  subsequent  fixes  can  use 
the  preceding  fix  as  their  approximate  point.  This 
method  will  work  fine  until  a  baseline  is  crossed,  at  which 
point  an  error  will  be  made.  To  prevent  this  error,  a  new 
approximate  location  must  be  specified  on  the  correct  side 
of  the  baseline.  In  most  cases,  this  is  taken  care  of 
automatically  if  the  approximate  point  is  extrapolated  from 
the  last  few  locations.  Still,  when  a  baseline  is  crossed, 
the  user  should  check  the  results  carefully. 

IV.  Conclusions 

The  methodology  described  here  has  been  implemented 
for  one  year  and  has  survived  an  initial  period  of  testing.  It 
is  part  of  the  NAVOCEANO  Hydrographic  Post  Time 
System  (HPTS).  Test  data  for  initially  debugging  the  pro¬ 
grams  is  presented  in  Appendix  II. 
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Subroutine  MPP 


(Given  an  approximate  position  and  a  set  of  navaid  observations,  compute  position.) 

Read  navaid  observations  and  depth. 

Begin  iterative  solution:  Repeat  until  positional  error  is  within  tolerance  or  the  number  of  iterations 
is  excessive. 

1  Do  2  I  =  1,  4  (once  for  each  navaid) 

If  aid  was  not  selected,  or  data  is  unavailable,  set  weight  to  zero,  go  to  2 
CASE  (R,  A,  H) 

RANGE 

Compute  the  range  and  azimuth  from  the  approximate  position  of  the  ship  to  this  navaid. 
Compute  B ,  f  and  W  matrix  elements. 

AZIMUTH 

Compute  the  range  and  azimuth  from  the  approximate  position  of  the  ship  to  this  navaid. 
Compute  B  matrix  elements. 

Compute  back  azimuth  (from  center  station  to  approximate  position  of  the  ship). 
Compute  f  and  W  matrix  elements. 

HYPERBOLIC. 

Compute  range  and  azimuth  from  approximate  position  of  ship  to  the  master  and  slave 
stations. 

Correct  ranges  for  long-range  propagation  over  salt  water. 

Compute  B ,  f  and  W  matrix  elements 
END  CASE 

2  END  DO 

(At  this  point,  all  matrices  are  filled.) 

Solve  Equation  (20)  for  position  correction  (dx,  dy). 

Apply  corrections  to  ship  position. 

If  magnitude  of  correction  is  within  tolerance,  RETURN. 

Else  if  number  of  iterations  is  not  excessive,  GO  TO  L 
Else  print  “NO  CONVERGENCE”. 

RETURN 


> 
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Appendix  II.  Test  data 


A.  Range -Azimuth 

The  following  data  were  used  to  test  combinations  of  range  and  azimuth  navaids.  The  correct  posi¬ 
tion  is  08°15T8.211'S,  116°57  11.205 "E. 

1.  Range  data 

Lanewidth  =  87  m 


Navaid 

Latitude 

Longitude 

Observation  (lanes) 

1 

8°14'23.0155  "S 

116°52'43.710  "E 

96.11 

2 

8°17'18.3105  "S 

116°55T7.110  "E 

58.40 

2.  Azimuth  data 

Navaid 

Station 

Latitude 

1 

Center 

8°14'23.125  "S 

Target 

8°16'38.080  "S 

2 

Center 

8°17'18.4515  "S 

Target 

8°16 ' 38.0805  "S 

Longitude 

Observation 

a0 

116°52 '43.937  "E 

317.370° 

0.01 

116°54'21.159"E 

116°55  17.151  "E 

97.479° 

0.01 

116°54'21.159  "E 

B.  Hyperbolic 

The  following  table  presents  data  for  5  positional  fixes. 

LORAN  A 
Clarke  1866  spheroid 


Latitude 

Longitude 

Coding  Delay 
(/xsec) 

Speed  of  Propagation 
(m//xsec) 

Master 

N41  14  56.330 

W69  58  31.4600 

Slave  1 

N35  14  25.9300 

W75  31  37.8300 

1000. 

299.692 

Slave  2 

N43  27  33.4500 

W65  28  16.3300 

1000. 

299.692 

Observations 
(time  differences,  fisec) 

Computed  Positions  (two  methods) 

TDl 

TD2 

GPTD  UNIVAC  1108 
O.S.  AA2NAA 

WITS 

Position  Difference 
(seconds) 

4400.00 

2800.00 

N35  24  03.7116 

N35  24  03.7112 

.0004 

W64  33  05.4840 

W64  33  05.4835 

.0005 

5800.00 

1900.00 

N39  56  47.1273 

N39  56  47.1270 

.0003 

W62  48  00.2974 

W62  48  00.2966 

.0008 

3900.00 

3300.00 

N35  37  49.0375 

N35  37  49.0372 

.0003 

W67  54  02.0548 

W67  54  02.0544 

.0004 

6000.00 

2800.00 

N40  23  02.8754 

N40  23  02.8752 

.0002 

W66  59  26.9214 

W66  59  26.9211 

.0003 

2400.00 

3800.00 

N35  26  49.4144 

N35  26  49.4137 

.0007 

W72  30  20.6275 

W72  30  20.6269 

.0006 

Source:  NAVOCEANO  Computer  Branch  (1972).  WITS  Program  Documentation. 

Naval  Oceanographic 

Office,  NSTL,  Mississippi. 
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